Cluster formation in compact stars: relativistic versus Skyrme nuclear models 
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dependent couplings) and non-relativistic Skyrme-type interactions. We mainly focus on the liquid- 
gas instabilities that occur at sub-saturation densities, leading to the decomposition of the homoge- 
neous matter into a clusterized phase. Such study is related to the description of neutron-star crust 
(at zero temperature) and of supernova dynamics (at finite temperature). 
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I. INTRODUCTION 

The knowledge of the equation of state (EOS) of nu- 
clear matter under exotic conditions is essential for our 
understanding of the nuclear force and for astrophys- 
ical applications. This implies high isospin asymme- 
tries, finite temperatures, and a wide density range (both 
for subsaturation and suprasaturation densities). The 
next generation of observational and experimental data 
is expected to bring new constraints in order to refine 
the theoretical models: for instance, the forthcoming 
radioactive-ion-beam facilities (such as FAIR@GSI and 
SPIRAL2@GANIL) will allow the investigation of the 
isospin degree of freedom in nuclear structure and dy- 
namics. 

The present work is dedicated to the predictions of dif- 
ferent effective nuclear models. It is mainly focused on 
the liquid-gas instabilities present in nuclear and stellar 
matter at sub-saturation density. These instabilities are 
directly related to the bulk EOS. They are used to ex- 
plain the multifragmentation phenomenon occurring in 
collisions around the Fermi energy [1]: in the spinodal 
decomposition scenario, fragment formation is induced 
by the fast development of spinodal instabilities in the 
low-density expanding matter formed just after the colli- 
sion [2] . Finite-size liquid-gas instabilities are also impor- 
tant for compact-star physics: matter non-homogeneities 
in the (hot) core of type II supernovae is expected to 
affect the dynamics of the explosion, and the crust of 
(cold) neutron stars contains a non-homogeneous phase 
commonly named pasta phase [3-5]. It should be no- 
ticed that the study of liquid-gas instabilities is comple- 
mentary to the equilibrium approaches which are also 
used to describe the clusterized stellar matter: namely 
(at very low density) the virial equation of state [6] , and 
(at higher density) the calculation of the pasta phases 
as the ground state shaped by the competition between 
Coulomb repulsion and surface tension. Although nu- 



clear equilibrium is expected to be reached in most stel- 
lar conditions, the spinodal-instability properties should 
help to understand the physics of compact stars in the 
following ways: (i) giving an estimation of clusterized- 
matter properties, such as cluster size and composition; 
(ii) showing the minimal region where the equilibrated 
matter must be formed of clusters; (iii) possibly, playing 
a direct role in cluster formation for specific situations 
involving very short time scales (as may happen during 
a supernova explosion). 

In this paper, we compare predictions from two kinds 
of models based on phenomenological density functionals: 
relativistic and non-relativistic. Both are commonly used 
to describe asymmetric nuclear matter, in the framework 
of exotic nuclei as well as compact stars. However, it is 
well-known that, although all give a quite good descrip- 
tion of stable nuclei (consistently with the constraints 
included in the fitting procedures), they present different 
behaviours as soon as exotic conditions are reached, es- 
pecially in the isovector channel. Our scope is to explore 
the impact of these different behaviors on quantities of in- 
terest for compact-star physics, such as the clusterization 
properties and the matter composition at /3-equilibrium. 

As a non-relativistic approach, we use the effective 
density dependent Skyrme-type interaction [7-9]. The 
simple form of the Skyrme functional makes it an at- 
tractive model for the description of both nuclei and 
compact-star matter. It was originally intended to de- 
scribe nuclear properties through the mass table, and 
the older parametrizations only include in their fits con- 
straints from magic-nucleus properties along the stability 
line. Trying to give a reliable description of exotic nuclei 
and stellar matter, the modern Skyrme parametrizations 
also include in their fitting procedures results from micro- 
scopic calculations of neutron-rich matter. The Skyrme- 
Lyon (SLy) forces for instance have been used in stud- 
ies of neutron-star crust [10, 11]. Such parametrizations 
are among the 27 forces which were not ruled out for 
unfit neutron-star properties in the extensive study by 
J.R. Stone et al. [12], where 87 Skyrme parametrizations 
were checked. 

In contrast with the Skyrme approach, the relativis- 
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tic nuclear models are, by construction, causal and can 
thus be applied to a wider region of the compact stars 
(as long as matter is supposed to be in an hadronic 
phase). Relativistic Mean-Field models (RMF) have 
been used to describe the EOS of compact stars [13-15], 
both cold and warm. The Density Dependent Hadronic 
models (DDH) are an alternative approach for the de- 
scription of nuclear matter and finite nuclei [16]. In 
DDH, the non-linear self-interactions of the mesons oc- 
curring in constant coupling models are substituted by 
density-dependent meson-nucleon coupling parameters, 
motivated by Dirac-Brucckncr calculations of nuclear 
matter. Such models are found to behave more closely 
to the non-relativistic ones. 

Relativistic and Skyrme approaches have been com- 
pared from the formal point of view in a recent work [17], 
where a low-density expansion of the RMF and DDH 
models have been used in order to directly compare the 
different density functionals. oth kinds of models have 
been used separately in several previous works for the 
study of spinodal instabilities in nuclear and compact- 
star matter, at zero and finite temperature: see for in- 
stance Refs. [18-22] for relativistic models, and Refs. [23- 
25] for Skyrme models. The same qualitative features 
are reproduced (general shape of the instability regions, 
isospin-distillation property of the phase separation). 
The scope of the present paper is then to have a direct 
look at the quantitative differences between relativistic 
and Skyrme-model predictions. We wish to investigate 
the extent to which the different temperature and isospin 
dependences of the nuclear EOS can affect the neutron- 
star properties and determine the sensitive features that 
have to be constrained. 

In section II we briefly review the relativistic and 
Skyrme models used in the present work. In section III, 
we present the Vlasov formalism that we use to address 
the dynamic instabilities, in both frameworks. Nuclear- 
matter properties are discussed for the different mod- 
els in section IV, where we present the nuclear EOS (in 
isoscalar and isovector channel) as well as the spinodal 
instabilities (in both thermodynamical and dynamical 
frameworks). Properties of stellar matter, including ho- 
mogeneous /3-equilibrium matter and instabilities against 
clusterization, are discussed in section V. In the last sec- 
tion we draw the main conclusions from our work. 



II. EFFECTIVE NUCLEAR MODELS 



A. Skyrme functional 

The local Skyrme interaction [8] allows to introduce an 
energy density H(r) so that the total energy for a system 
of nucleons in a Slater determinant | ip > reads : 



W\H\- 



H(r)d 3 r , 



(1) 



where H(r) is the Skyrme energy-density functional. 

In the case of homogeneous, spin-saturated matter 
with no Coulomb interaction, the Skyrme energy-density 
functional [28] reduces to four terms: 



H = JC + Ho + H 3 + Tt c fi 



(2) 



where the label b (bulk) is used to mark the thermody- 
namic framework. In this expression, JC is the kinetic- 
energy term, Ho a density-independent two-body term, 
Hz a density-dependent term, and H e s a momentum- 
dependent term: 



JC 
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H 3 
H c s 



T 

2m 

Cop 2 + Dopl 
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(3) 

(4) 
(5) 
(6) 



We have introduced the isoscalar and isovector particle 
densities, p and p 3 , as well as kinetic densities, t and t 3 : 



P = Pn+ Pp 
P3 = Pn - Pp 



T = T n + T p 
T 3 = T n - T p 



(7) 



where, denoting i the third component of the isospin (n 
for neutrons and p for protons) , the kinetic densities are 
defined by n — (k )i. The coefficients C and D, asso- 
ciated respectively with the symmetry and asymmetry 
contributions, are linear combinations of the traditional 
Skyrme parameters: 



Co — 3io/8 

D = -t (2x + l)/8 

C 3 = is/16 

£> 3 = -t 3 (2a;3 + l)/48 

Cos = [3t 1 +t 2 (4 a ; 2 + 5)]/16 

D eS = [t 2 (2x 2 + l)-t 1 (2x 1 + l)}/16 



(8) 



In the sequel we will give a short presentation of the 
models discussed in the present paper. All expressions 
are given in units h = c = 1. We will consider first den- 
sity functionals based on Skyrme forces, then the RMF 
and DDH models. The nuclear matter saturation prop- 
erties obtained with all models used in the present work 
are reported in Table I. 



In the mean-field approach, the individual particle level 
is derived from this functional. For each particle species, 
it is given by: 



dH b dH b 



dpi 



+ 



dn 



k 



= mi + U l + 



2m* 



-kf. 



(9) 
(10) 



3 



TABLE I: Nuclear matter properties of the Skyrme and relativistic models used in the present work. 



Model 


B/A 


Po 


K 


ra* /m 


a a (po) 


L(po) 


^Csym(po) 




(MeV) (mi" 3 ) 


(MeV) 




(MeV) (MeV) 


(MeV) 


SIII [26] 


15.9 


0.145 


356 


0.76 


28.2 


9.9 


-394 


SGII [27] 


15.6 


0.159 


215 


0.79 


26.9 


37.6 


-146 


SLy230a [28] 


16.0 


0.16 


230 


0.70 


32.0 


44.3 


-98 


NRAPR [29] 


15.9 


0.16 


226 


0.70 


32.8 


59.6 


-123 


LNS [30] 


15.3 


0.175 


211 


0.83 


33.4 


61.5 


-127 


NL3 [31] 


16.3 


0.148 


272 


0.60 


37.4 


118.3 


101 


NL<5 [32] 


16.0 


0.160 


240 


0.75 


30.5 


102.7 


127 


TW [33] 


16.3 


0.153 


240 


0.56 


32.0 


55.3 


-125 


DD-ME2 [34] 


16.1 


0.152 


251 


0.57 


32.3 


51.7 


-88 


DDH<5 [35] 


16.3 


0.153 


240 


0.56 


25.1 


48.6 


81 



where we have included the nuclear mass energy m. The 
kinetic energy is expressed in the non-relativistic limit, 
in terms of an effective mass m* defined by: 



2m* 2m 



1 + dH eS 



(11) 



The chemical potentials fii are such that the Fermi- 
Dirac occupation number is: 



n,i{k) 



1 + e 



(12) 



we can also define a chemical potential \i\ = [n + mi 
including the mass energy, such that: 



rii(k) = 



1 + e 



(13) 



In this work, we will use conventional and modern 
Skyrme interactions. The earlier parametrizations, such 
as SIII [26] and SGII [27], have been established by fit- 
ting the properties of stable nuclei (such as radii and 
ground-state energy). They are thus in principle poorly 
adapted to a description of neutron-rich matter. It is 
indeed found that SIII presents an irrealistic behavior 
in the isovector channel; furthermore it has a too high 
incompressibility at saturation. As a result, this inter- 
action will present an atypical behavior all through the 
following study. SGII, however, for which spin proper- 
ties have also been used as constraints, presents a more 
reasonable evolution in the isovector channel. In par- 
ticular, it has been shown to reproduce isospin effects 
in giant dipole resonances [27]. Among the modern 
Skyrme-type parametrizations, we have chosen to use 
one of the Skyrme- Lyon forces (SLy230a [28]), as well 
as the NRAPR [29] and LNS [30] parametrizations. All 
these recent forces include in their fitting procedure re- 
sults from microscopic calculations. SLy230a uses the 
pure-neutron matter equation of state UV14+UVII by 
R.B. Wiringa et al [36]. NRAPR (Non-Rclativistic APR) 
stands for the Skyrme interaction parameters obtained 



from a fitting to the APR equation of state (Akmal- 
Pandharipande-Ravenhall, Ref. [37]). LNS is based on 
Brueckner-Hartree-Fock calculations of infinite nuclear 
matter at different values of isospin asymmetry. Such 
constrains from microscopic calculations are intended to 
control the behavior of the resulting effective force far 
from saturation and up to high isospin asymmetry. 



B. Relativistic approaches 

In the present paper, we will consider two kinds of 
relativistic effective approaches: RMF models, which 
have constant coupling parameters described by the La- 
grangian density of non-linear Walecka models (NLWM) , 
and DDH models with density-dependent coupling pa- 
rameters. In each case, we consider models including or 
not the (5-meson, which have been introduced to include 
in the isovector channel the same symmetry existing al- 
ready in the isoscalar channel with the meson pair (ct, uj) 
responsible for saturation in RMF models [32]. The pres- 
ence of the 5-meson softens the symmetry energy at sub- 
saturation densities and hardens it above saturation den- 
sity. The RMF parametrizations we use are NL3 [31] and 
NL(5 [32]; the DDH ones are TW [33], DD-ME2 [34] and 
DDH<5 [35]. Only NL(5 and DDHJ include the 6-meson. 

The relativistic approach is based on a lagrangian den- 
sity given by: 

C = Yl Ll +C*+Cu+C p + C s . (14) 

i—p.n 



The nucleon Lagrangians read: 



with 



A = V>i bpiD* - M*] V>, 



iD» = id* -r„V" - 

2 

M* = m-T s (j)-T 5 T- 5 , 



(15) 

(16) 
(17) 
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where f is the isospin operator. We use the vector symbol 
to designate a vector in isospin space. 

The isoscalar part is associated with the scalar sigma 
(a) field cj> and the vector omega (ui) field V p , while the 
isospin dependence comes from the isovector-scalar delta 
(6) field 5 l and the isovector- vector rho (p) field 6^ (where 
p is a space-time index and i an isospin-direction index). 
The associated Lagrangians are: 



r 

Co 



1 



-m 2 s 4> 2 )-^-^ 



2p 
5 C 

1 



where = d^V v - d^V^, B pv = d^K - d v b p - T p (b p x 
b v ), and Tj and rrij are respectively the coupling param- 
eters of the mesons j — s,v,5,p with the nucleons and 
their masses. The self-interacting terms for the cr-meson 
are included only for the NL3 and NL<5 parametrizations, 
k and A denoting the corresponding coupling constants. 

The density-dependent coupling parameters r s , r„ and 
T p , are adjusted in order to reproduce some of the nuclear 
matter bulk properties, using the following parametriza- 
tion: 



r i(/>) = ^i{Psat)fi{x) , i = s,v 



with 



fi(x) = a, 



l + b l {x + d l )' 2 



(18) 



(19) 



' \ + Cl (x + d t Y ' 

where x — pj p sa t and 

T p (p) = T p { Psat ) exp[-a p (z - 1)] . (20) 

The values of the parameters mi, r,, aj, bi, Ci and di, 
i = s,v,p for TW and DD-ME2 are respectively given 
in [33] and [34] and for DDH<5 in [18, 35]. In this last case 
the parametrization for the 6 and p coupling parameters 
is also given by (18) with 

fi(x) = a,i exp[-6i(a; - 1)] - a(x - di) , i = p, 6. 

The Ti coupling parameters are replaced by the gi cou- 
pling constants in the NL3 and NL<5 models. 



Skyrme and relativistic models. For simplicity, we con- 
sider here nuclear matter where the proton electric charge 
is neutralized in average by a uniform background. We 
thus neglect the electron degree of freedom existing in 
star matter, which is model-independent and has only a 
perturbative effect. 



A. Brief review 

To describe the time evolution of the nuclear system, 
we introduce the one-body phase-space distribution func- 
tion in isospin space: /(r, k, t) = diag(/ p ,/„), and the 
corresponding one-body Hamiltonian h = diag (h p , h n ). 
The time evolution of the distribution function is de- 
scribed by the Vlasov equation: 



dt 



+ {fi,hi} = 0, i=P,n, 



(21) 



where {,} denotes the Poisson brackets. At zero tem- 
perature, the state which minimises the energy of asym- 
metric nuclear matter is characterised by the Fermi mo- 
menta kpi, i — p,n, and is described by the distribu- 
tion function / (k) = diag [Q{k 2 Fp -k 2 ), Q(k 2 Fn ~k 2 )]. 
In order to describe small oscillations around the equi- 
librium state, we take for the distribution functions 
fi = foi + Sfi an d introduce a generating function [38] 
S(r, k, t) = diag(S'p, S n ) defined in isospin space such 
that 5fi = {Si,foi}. In terms of the generating function, 
the linearised Vlasov equations for 5fi are equivalent to 
the following time evolution equations: 



dS- 

+ {Si,h 0i } = {5hi) F , 



(22) 



where (5hi)p is the mean- field variation at Fermi level, 
which depends on the considered nuclear model. 

We will consider the longitudinal fluctuations such 
that: 



( Si; Sp l ; Shi^j - ( S U! i(x); 5pu>,i\ ) e l(q ' r ut) 

; (23) 

where x = cos(k, q). The longitudinal normal modes are 
obtained substituting the ansatz (23) in the linearised 
equations of motion. The dispersion relation takes the 
form: 



1 + FPPL P FP n L p 



A, 



, , ^ = 0, (24) 
F n P L n 1 + F nn L n \A un 



III. THE VLASOV FORMALISM 

In this paper, we study the dynamic spinodal insta- 
bilities as unstable density fluctuation modes obtained 
in the Vlasov framework. The present section gives a 
short review of the Vlasov formalism already introduced 
in Refs. [20, 38, 39] and give the resulting expressions for 



where Li = L(si) = 2 — Sjln[(sj + l)/(sj — 1)] is the 
Lindhard function, with Sj = uj/(qvFi) in terms of the 
Fermi velocity VFi — deFi/dkFi- 

The amplitudes A^i = J^x S u i(x) dx are related to 
the transition densities by 

_ uNpj A 

"Pi — — A*n; 

2 Si 



■5 



where Noi is the density of states at the Fermi sur- 
face. With all models, the coefficients F y appearing in 
Eq. (24) can be expressed in terms of two quantities re- 
lated to the nuclear residual interaction, U-P and uf^ 
defined by: 

{Shi) F = ^[uW+xUg^Spj. (25) 



In this expression, we separate the x-dependent contri- 

/(2) 



bution of the residual interaction {xU\j ) from the x- 

independent one (U^p). On the other hand, Eq. (22) 
provides the following relation: 

{5hi) F = -iwSi(x) [1 - x/si] . (26) 



From Eqs. (25) and (26), we obtain the set of Vlasov 
equations: 



N . 



■6 Pi 



(2) 



5 PJ = , (27) 



which is equivalent to Eq. (24) with the following identi- 
fication: 



N 0j s, f 



(28) 



At zero temperature we have Srj/Spj — kp^, and taking 
the value at Fermi level we get: 



{Sh.% 



E 



d 2 n d 2 n 2 2 ' 

d Pj dpi + a^dp^ + kFt \ 



S Pj (32) 



The surface term arises from the density-gradient de- 
pendence in the Skyrme Hamiltonian density: 

n v = cJ„(v P „) 2 + cJ p (v Pp ) 2 + 2c„ v p v Pn v Pp ) 

= CE(Vp) 2 + C£(Vp 3 ) 2 , (33) 

where the coefficients C v are combinations of the usual 
Skyrme parameters (given in Rcf. [40]), independent of 
neutron and proton densities. A transferred momentum 
q then induces the nuclear mean-field variation: 



ShJ = 2q 2 ^C7 j 5p j . 



(34) 



Let us finally consider the Coulomb mean-field varia- 
tion. In the non-relativistic limit, for nucleons with ef- 
fective mass to* we have: 



w)f - E 



^TTCiCj 1 — xv* Fi uj/q 



l-oj 2 /q 2 



Spj , (35) 



The model-dependence of the Vlasov equations is then where v* Fi = k Fi /m*; e„ = 0; e p = e = q e /y/4^. Only 



(1 2) 

contained in the cocfficcnts W- ' that we have intro- 
duced: next, more details will be given for both Skyrme 
and relativistic models. 

Let us note that the instabilities of the system are de- 
termined from the imaginary frequencies which satisfy 
the dispersion relation [39] . The finite-size instability re- 
gion is the envelope of all the dynamical spinodals corre- 
sponding to different values of the transfered momentum 
q. For simplicity, through out this work we will identify 
this domain with the dynamical spinodal for q = 80 MeV, 
which is a very good approximation for all the models un- 
der study. 



the Coulomb term brings a x-dependence in Shi 
From the above expressions, we identify: 



*3 



d 2 n 



d 2 n 



dp J dp l + drjdp, ^ 



Fj 



4ireiej 



V& = 



q 2 1 — Ct> 2 /g s 
47reje 3 v* Fi uj/q 
q 2 1 — uj 2 /q 2 



(36) 
(37) 



B. Vlasov approach with Skyrme models 

We identify three contributions to the mean-field vari- 
ation Sh^. the bulk (6), surface (V) and Coulomb (c) 
terms, such that 



Shi = 6h b i + Shf + Sh c i 
The bulk term is: 

k 2 



5h\ = 5 

= E 



Ui 



2m* 
d 2 H 



(29) 



(30) 



0-H .St. 



dp 3 dpi drjdpi Sp 



k 2 ) 



5 Pj . (31) 



C. Vlasov approach with relativistic models 



For relativistic models the one-body Hamiltonian is 
written in terms of the meson fields 

hi = \J(k-Vi) 2 +mf + Va, i = p,n, 

where m* = m — Y s <p — T^gS^ denotes the effective mass 
of nucleon i and 

v 0l = r v v + ^ nb + eAo^p + , 



V i =T v V+'^-T i b + eA^i, 
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with n = 1 (— 1) for protons (neutrons). The contribu- 
tion of the rearrangement term, due to the density de- 
pendence of the coupling parameters is given by 

R dr v dv p b dr s dr s 

The variations of the one-body Hamiltonian which en- 
ter the linearised Vlasov equations are: 

Shi = 6(m* - m) 1 ^ + SV i - , (38) 

with 

h 0i = yJv + mf + V™ = eio+V® 

and 5{m*-m) = —(T s 8(f) + 8T s (po+n Tg SS 3 + n S3 8Tg). 
The x-dependence of Shi is present through the contri- 
bution of the spatial components of the vector fields: the 
to and p mesons and the electromagnetic field. 

Using the linearised equations of the fields, we express 
the field variations in terms of the proton and neutron 
particle densities and scalar densities [39] and reduce Shi 
to an expression similar to (25). The coefficients i™ have 
been defined in [20, 22]. The nuclear-energy dependence 
on the transfered momentum involves the different meson 
masses, and is more complex than the Skyrme quadratic 
expression. This point will be discussed in more details 
in Sec. IV D. 



IV. NUCLEAR MATTER PROPERTIES 

In the present section we will compare the nuclear mat- 
ter properties predicted by all the models under study. 
We shall consider first the isoscalar properties of the 
EOS, then the isovector ones. The spinodal instabilities 
will also be analysed: we will address the thermodynamic 
instability region and direction of phase separation, and 
finally the clusterization properties within the Vlasov ap- 
proach. 

A. Symmetric Nuclear Matter 

In Fig. 1 we show some bulk isoscalar properties of 
nuclear matter as a function of the baryonic density, 
namely a) the energy per nucleon, b) the pressure P = 
p 2 d(E/A)/dp, and c) the incompressibility K — 9dP/dp. 
Curves are shown for Skyrme interactions (left) and rela- 
tivistic models (right). Skyrme interactions show similar 
behaviours between them, except for SIII (stiffer) and 
LNS (softer). 

Globally, the relativistic models present slightly higher 
binding energies, lower saturation densities, higher in- 
compressibilities (disregarding SIII) and lower effective 
masses; however, we can note that NL<5 presents isoscalar 
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FIG. 1: (Color online) Bulk isoscalar properties of nu- 
clear matter, as a function of the baryonic density for non- 
relativistic (left) and relativistic (right) models. From top to 
bottom: energy per nucleon, pressure and incompressibility. 

saturation properties quite similar to the Skyrme ones. 
Let us remind that the effective mass has a different 
meaning in each framework, as was already stressed in 
other works [17, 41]: in relativistic models, it includes 
the contribution of the nucleon scalar self-energy, while 
for the Skyrme interactions it reflects the momentum de- 
pendence of the single-particle energy. 

B. Asymmetric Nuclear Matter 

We will now discuss properties of asymmetric nuclear 
matter (ANM). In Fig. 2 we show some properties re- 
lated to the isovector channel of the respective models. 
As expected, larger differences are observed in this chan- 
nel; not only between relativistic and Skyrme models, 
but also between different parametrizations inside each 
framework. We plot the symmetry energy 

_ 1 d 2 {E/A) _ 1 d 2 (E/V) 
s 2 dy 2 2p dy 2 ' 
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FIG. 2: (Color online) Bulk isovector properties of nu- 
clear matter, as a function of the baryonic density for non- 
relativistic (left) and relativistic (right) models. From top 
to bottom: the symmetry energy and its derivatives with re- 
spect to the density, namely the slope parameter L and the 
symmetry incompressibility K syln . 

with y = (p n — p p ) I (p n + p p ) = pzl P- Note that this defi- 
nition can be generalised at finite temperature: denoting 
T the free-energy density, it becomes 

_ 1_ d 2 T 

fa 2p dy 2 

Figure 2 also represents quantities related to the first and 
second density derivatives of the symmetry energy (re- 
spectively denoted by a' s = da s /dp and a" = d 2 a s /dp 2 ), 
according to expressions of common use in the litera- 
ture [42]: the slope parameter 

L = 3p a' s 

related to the symmetry pressure at saturation, and sym- 
metry incompressibility 



of the symmetry energy in the presented density range. 
As expected, the older SIII parametrization presents 
atypic features; it even predicts an isospin instability at 
p = 0.325 fm~ 3 , as it can be seen in Fig. 2 top. SGII 
follows an intermediate behaviour. It is interesting to 
see that NRAPR and LNS, despite different symmetric- 
matter EOS, almost coincide in the isovector channel. 
With L values at saturation of the order of 60 MeV, these 
two parametrizations get close to the L range estimated 
from the most recent experimental constraints (isospin 
diffusion and isoscaling data): L — 88 ± 25 MeV [43- 
45], while the other three Skyrme paraterizations have 
too small symmetry-energy slopes at saturation. The 
parametrization SLy230a differs from NRAPR and LNS 
by its sharp increase of the symmetry energy at suprasat- 
uration densities. In the following we will see how these 
quantities influence the predictions of the different mod- 
els for neutron rich matter. 

Among the relativistic parametrizations used which do 
not include the 5-meson, the larger differences occur be- 
tween NL3 and the models with density dependent cou- 
plings: TW and DD-ME2, which have very similar be- 
haviours in the isovector channel. In particular, NL3 
has a very hard symmetry energy which increases almost 
linearly with the density. The inclusion of the <5-meson 
significantly reduces the symmetry energy at p < po, but 
it is also associated with a sharp increase of a s at higher 
densities. Considering the L value at saturation, we see 
a clear separation (roughly a factor of 2) between RMF 
and DDH models, situated on each border of the interval 
of experimental constraints cited above. 

It is interesting to see that TW and DD-ME2 (DDH) 
behave like NRAPR and LNS (Skyrme) for all isovec- 
tor properties. However, all other parametrizations show 
large differences affecting the three quantities a s , L and 
K sym . The general trend is that relativistic models have 
a stiffcr symmetry energy, as well as a larger symme- 
try incompressibility. For the presented results, several 
relativistic models have a region of positive K sym , while 
SLy230a is the only Skyrme parametrization to present 
such feature. In the following, we will try to investi- 
gate to what extent the /3-equilibrium and clusterization 
properties are affected by these differences. 

Let us finally investigate the validity of the parabolic 
approximation of the isovector EOS, which is model de- 
pendent. In this approximation, we have a direct link 
between the symmetry energy and the isovector chemi- 
cal potential /13 = p n — p p , which determines the matter 
composition at /3-equilibrium. Indeed, the parabolic ex- 
pression of the free-energy density is: 

■F~jr paia = T s + pa s y 2 (39) 

where T s = !F(p, 0) is the free energy density of symmet- 
ric matter. The corresponding isovector chemical poten- 
tial is then proportional to y, as: 



Among the Skyrme forces, the modern parametriza- 
tions (SLy230a, NRAPR and LNS) show similar values 



(40) 
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FIG. 3: (Color online) Isovector chemical potential p- A = fi n — p p as a function of the asymmetry y = ps/p for two values of 
the total density (p = 0.05 and 0.3 fm -3 ). Top: exact value of p-j,. Bottom: the ratio ps/p^^, where p^"" — 4a s pzj p- 



The parabolic approximation is exact in the limit of small 
asymmetry, and actually gives very good predictions for 
T until y = 1. However, more significant differences may 
be obtained for the pa-derivative leading to /13. This be- 
haviour is checked on in Fig. 3, representing both the ex- 
act value /13 = 2dJ 7 /dp3 and the ratio p3/pg ara . For this, 
we have fixed two values of the baryonic density, p = 0.05 
and 0.3 fm~ 3 . At low densities we mostly confirm the va- 
lidity of the parabolic approach; more significant differ- 
ences are observed at high density. The dominant trend 
is to have ^3 > p^ aia , due to the kinetic contribution to 
the symmetry energy; only the SLy230a and DDH<5 show 
the opposite behaviour, at high density. 



C. Thermodynamical spinodal instability 

The liquid-gas phase transition is a well-known fea- 
ture of the nuclear-matter EOS. It corresponds to the 
presence of an abnormal (negative) curvature of the free- 
energy density T as a function of (p„, p p ), or equivalently 
(p, P3). The thermodynamical spinodal instability cor- 
responds to the region where the homogeneous matter 
is locally unstable against the separation in two infinite 
homogeneous phases, meaning that the surface T(p,pz) 
presents a local negative curvature. This bulk property 
of nuclear matter is at the origin of the dynamic insta- 
bilities leading to matter clusterization. We will consider 
next the thermodynamical spinodal properties. 



Thermodynamic spinodal region 



The spinodal contour is defined by the cancellation of 
the determinant of the free-energy curvature matrix: 



C 




dpidpj 



(41) 
(42) 



where p\ = p. Inside the spinodal region, the lower eigen- 
value C< of this matrix is negative. 

The different spinodal contours are shown on Fig. 4, 
using two different representations: the density plane 
(p n , p p ) and the mixed plane (p, ^3), where isovector dif- 
ferences appear more clearly. Although the overall fea- 
tures are similar, the trend is that relativistic models 
predict a smaller instability region, both in isoscalar and 
isovector directions. The isoscalar extension of the spin- 
odal is measured by the density p s , corresponding to the 
upper spinodal border for symmetric matter. The dif- 
ferent p s values are reported in Table II: we can verify 
that they are correlated with the po values. As for the 
isovector behavior of the spinodal contour, we see that it 
reaches very high asymmetries with all models. We can 
however compare the different extensions obtained in the 
P3 direction. They are found to reflect the subsaturation- 
density behavior of the symmetry energy: indeed, as dis- 
cussed above, in this density range the isovector chemical 
potential can be well- approximated by ps ~ 4et s ?/. For 
instance, NL5 has the smallest a s values at low density, 
and therefore presents the narrowest spinodal contour in 
the P3 direction. We remind that the inclusion of the 5- 
meson leads to a reduction of a s , which is observed both 
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TABLE II: Characterization of the spinodal shape by the con- 
tour concavity C B , depending on the symmetry energy and its 
density derivatives taken at the upper spinodal border of sym- 
metric matter (p — p s ). 





Po 


Ps 


o s 


p s a' B 


2 // 

P s a s 






(frn" 3 ) 


(fm" 3 ) (MeV) (MeV) (MeV) (MeV fm 3 ) 


SIII 


0.145 


0.098 


24.73 


12.06 


-20.96 


-175.72 


SGII 


0.159 


0.100 


21.06 


12.41 


-9.33 


17.15 


SLy230a 


0.160 


0.102 


25.73 


13.21 


-9.70 


61.60 


NRAPR 


0.161 


0.103 


24.60 


16.50 


-8.24 


51.22 


LNS 


0.175 


0.111 


24.89 


16.92 


-8.28 


46.27 


NL3 


0.148 


0.096 


24.04 


23.40 


2.86 


85.98 


NLS 


0.160 


0.102 


19.06 


18.93 


3.87 


81.33 


TW 


0.153 


0.096 


24.71 


15.95 


-9.91 


29.10 


DDME2 


0.152 


0.099 


25.57 


15.00 


-10.92 


29.98 


DDH<5 


0.153 


0.096 


19.80 


10.38 


-6.03 


80.25 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



P 3 /P P 3 /P 




-1 -0.5 0.5 1 -1 -0.5 0.5 1 



P3'P P 3 'P 



with NL<5 and DDHJ. The ^-extension of the spinodal 
contour is a feature especially relevant in the astrophys- 
ical context, for the comparison between the instability 
region and the constraint of /^-equilibrium: this point will 
be addressed in more details in section V. 

Let us now consider the shape of the spinodal contour: 
the differences we observe can be caracterized by the con- 
vexity of the upper border. Therefore, we introduced 
the quantity C s (hereafter called contour concavity), de- 
fined as the convexity of the spinodal contour at point 



FIG. 5: (Color online) Thermodynamic instability direction 
as a function of the asymmetry, for two fixed total densi- 
ties. Results are given for Skyrme (left) and relativistic (right) 
models. 



(p = p s ,y = 0) [40]: 

9 2 C<(p s ,0) 



dy 2 

[pX + 2pX(l-pX/<>] ■ (43) 
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FIG. 4: (Color online) Thermodynamic spinodal of infi- 
nite nuclear matter for non-relativistic (left) and relativis- 
tic (right) models. The top spinodals are represented on the 
{p P ,Pn) plane and the bottom ones on the (p3,p) plane. 



If C s is positive (negative), for a small asymmetry y the 
point (p s ,y) is outside (inside) the spinodal, meaning a 
concave (convex) contour. Equation (43) gives the rela- 
tion between C s and the density behavior of the symme- 
try energy, involving a s , a' s and a"(p). None of the terms 
constituting this expression dominates, as can be appre- 
ciated in Table II: the symmetry energy, but also its first 
and second derivatives come into play to determine the 
contour concavity. Positive C s values are obtained with 
all the present models except SIII, whose convex shape 
is due to the large negative values of a"(p s ). Concerning 
the relativistic models, the models with constant cou- 
plings (NL3 and NL<5) have larger a' s and positive a" (p), 
giving rise to larger C s values. We can note that the 
contour concavity is a relevant property for the study 
of non-homogeneities in star matter, since it determines 
the sensitivity of the upper spinodal border to the specific 
composition that will be imposed by the /3-equilibrium. 



2. Thermodynamic instability direction 

The thermodynamic instability direction is the direc- 
tion of minimal free-energy curvature, given by the eigen- 
vector of matrix (41) associated with C<. It is related 
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FIG. 6: Thermodynamic instability direction as a function of 
the density, at proton fraction Y p — 0.3, for non-relativistic 
(left) and relativistic (right) models. 



to the phenomenon of isospin distillation, which usually 
leads to the formation of a dense phase more symmetric 
than the dilute one. We express this direction as the ratio 
5ps /Sp, giving the deviation with respect to the isoscalar 
direction. The eigen-vector (Sp, Sp^) < satisfies: 



Spa 
Sp 



C< — Tn 



13 



(44) 



This ratio is zero in the case of symmetric matter, where 
the instability direction is purely isoscalar. For ex- 
tremal asymmetry y = ±1, it obeys the limit conditions 
Sps/Sp — ±1, which constrains the behavior oiSp^/Splij) 
at high asymmetry. However, for moderate values of 
asymmetry, the evolution of the instability direction is 
nearly linear with respect to y, as illustrated in Fig. 2. 
In this region, the isospin-distillation properties of the 
different effective forces can be characterized by a num- 
ber 5 such that 



Sp3 



h + 0{y 3 ) 

a s - pa' s 
a s - pT></2 



(45) 
(46) 



On Fig. 6, we show the density evolution of Sp^/Sp 
at a fixed proton fraction Y p — 0.3, corresponding to 
the asymmetry y = p 3 /p = 0.4. Note that this re- 
flects the d values according to the good approximation 
d ~ (Sp3/Sp)/y. For all models, Sp^/Sp is lower than 
y = 0.4: this is the normal distillation effect. It is seen 
that the relativistic models with constant coupling, in- 
dependently of containing or not the (5-meson, predict a 
much larger distillation effect (smaller ratio Sp^/Sp): this 
is due to the quasi- linear behavior of a s (p) (see Fig. 5), 
which leads to low values of d. In contrast, the DDH 
models behave like the Skyrme forces: both show a re- 
duction of the distillation effect with density, while NL3 
and NL5 present the opposite behavior (which was also 
noticed in [19, 46]). 



D. Dynamical spinodal instabilities 



The bulk liquid-gas instability properties we have dis- 
cussed do not manifest themselves directly in nuclear 
multifragmentation and compact-star matter, due to the 
role of the Coulomb interaction and the surface tension. 
However, they induce instabilities against finite-size den- 
sity fluctuations, leading to the decomposition of the ho- 
mogeneous matter into a clusterized medium [47]. We 
now use the formalism presented in section III to study 
the Vlasov unstable modes, considering plane-wave den- 
sity fluctuations of wave-number q. The dispersion re- 
lation is defined by Eq. (24), with uj = i/r; t is the 
time constant which characterizes the initial growth of 
the density fluctuation. 

On Fig. 7, we compare the unstable modes obtained 
within the different models for fixed average densities 
(p = 0.05 fm~ 3 , Y p = 0.3), as a function of the wave num- 
ber. The top figures represent the growth rate |u/| = 1/t, 
and the bottom ones the direction of the mode in the 
density plane given by Sp^/Sp. 

Considering the top part of Fig. 7, we first see that 
relativistic models are usually characterized by a re- 
duced instability. To the noticeable exception of NL5, 
both the growth rate and the upper border of the unsta- 
ble q interval are smaller within the relativistic models. 
The bottom part of Fig. 7, gives the phase-separation 
direction Sp^/Sp associated with the dynamical modes. 
All curves decrease with q, due to the Coulomb effect: 
at low q, the strong Coulomb contribution quenches 
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FIG. 7: (Color online) Finite-size instabilities: dispersion re- 
lation for the unstable modes at p = 0.05 fm~ 3 , Y p — 0.3, 
T — 0, for Skyrme (left) and relativistic models (right). Top: 
modulus of the imaginary frequency \uj\ as a function of the 
wave number. Bottom: direction of the unstable mode Sps/Sp 
in the density plane. The dotted lines show pz/ p = 0.4 corre- 
sponding to the proton fraction considered. 
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FIG. 8: Most unstable modes for np matter with proton frac- 
tion Y p = 0.3, for Skyrme parametrizations (left) and rela- 
tivistic models (right). Top: growth rates. Bottom: associ- 
ated cluster size. 



the proton-density fluctuation, imposing large values of 
5 Pi 1 8 p. This leads to the so called anti- distillation effect, 
corresponding to Sp^/Sp > ps/ p, namely a dense phase 
more neutron- rich than the homogeneous matter: it is 
obtained here below q ~ 30 MeV/c. For higher q val- 
ues, we recover the normal distillation effect with all the 
models. Comparing the present dynamic results with 
the bulk instability direction (Fig. 6), we see that the 
behavior of the relativistic models deserves a comment. 
Although NL3 remains the model with the strongest dis- 
tillation effects, the hierarchy of the other curves is widely 
rearranged when dynamical instabilities are considered: 
the weaker distillation effect is now obtained with models 
including the 5 meson, both NLiJ and DDH5. 

Let us now study the most unstable mode, namely the 
mode of largest growth rate, which drives the system 
to the non-homogeneous phase. The associated wave- 
number characterizes the size of the primary clusters 
formed in spinodal decomposition, which we can define 
as the half wavelength of the fastest amplified mode. In 
Fig. 8 we show the growth rates and associated cluster 
size of the most unstable modes. We show results for 
Y p = 0.3, a proton fraction close to that of /3-equilibrium 
matter with neutrino trapping and to the asymmetry val- 
ues that could be involved in future multifragmentation 
experiments with radioactive beams (for instance 132 Sn 
has Y p = 0.379). Consistently with the observations of 
Fig. 7, as a general rule the Skyrme parametrizations 
predict larger growth-rates and smaller clusters than the 
relativistic models. The hierarchy between the differ- 
ent parametrizations is also essentially conserved: among 
Skyrme forces, SLy230a gives the largest clusters and 



LNS the smallest ones; among the relativistic models, 
NL<5 predicts particularly small sizes and DD-ME2 gives 
the largest clusters. Furthermore, we can notice features 
appearing with the density evolution. Firstly, going to 
lower densities, the Skyrme cluster sizes decrease more 
neatly than the relativistic ones. Thus, the minimal clus- 
ter sizes obtained in each framework are of ~ 2.5 fm 
(Skyrme) and ~ 4 fm (relativistic models). Secondly, 
larger sizes are reached near the border of the unsta- 
ble region (~ 8 fm with Skyrme, and beyond 10 fm 
with the relativistic models). These features accentuate 
the trend according which the relativistic models predict 
larger clusters. 

The different g-dependences originate in the finite- 
range part of the nuclear force, which was introduced 
in Sec. III. Since the direction of the density fluctuations 
is essentially isoscalar, we can characterize the energy 
cost of the density gradient by the quantity C^q 2 (for 
the nuclear contribution). This quadratic expression is 
exact for Skyrme models, but the q-dependence is more 
complex for relativistic models: performing a Taylor ex- 
pansion in powers of fi = q 2 /mf (m^ denoting the meson 
masses), we obtain a density-dependent CY\ coefficient. 
The values of are listed in Table III: for relativistic 
models, they are given at p — 0.05 fm -3 , together with 
the /, values. 

The correspondence between CYi and the maximal un- 
stable momentum q max is shown on Fig. 9 for symmet- 
ric matter, as a function of the density. Neglecting the 
Coulomb interaction (whose contribution is vanishing for 
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FIG. 9: (Color online) CYi and ~ maj t for Skyrme (left) and rel- 
ativistic (right) models, for symmetric matter. For relativistic 
models we compare g ma x,quad (thick lines) obtained from the 
quadratic expansion with the exact value q max (thin lines). 
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TABLE III: Dependence of the nuclear energy on the transfered momentum characterized by the CYi coefficients. These 
coefficients are constant for Skyrme models. For relativistic models, they are given at p = 0.05 fm -3 , together with the 
corresponding parameters: / 4 = g 2 /m 2 for i = a, u>, S, and f p = g 2 /(4m p ) 2 . 



Skyrme 




Relativistic 


c7i 


fs fv fp fs 




(MeV.fm 5 ) 




(MeV.fm 5 ) (fm 2 ) (fm 2 ) (fm 2 ) (fm 2 ) 


SIII 


63.0 


NL3 


99.2 


15.73 10.53 1.34 


SGII 


54.8 


NLS 


43.0 


10.33 5.42 3.15 2.5 


SLy230a 


77.7 


TW 


115.9 


18.97 14.64 1.79 


NRAPR 


64.1 


DD-ME2 


107.7 


18.50 13.99 1.94 


LNS 


43.8 


DDH<5 


115.9 


18.97 14.64 4.16 2.96 



the q range of interest), the quadratic expression of the 
(^-dependence leads to the following relation: 

n 2 = ^ (A7\ 

ymax.quad 2(7^ ' 

For the relativistic models, both g ma x,quad and the exact 
9max values are shown on the figure: the exact values are 
smaller than the values calculated in the quadratic ap- 
proximation. The different CYi values are seen to explain 
the different cluster sizes obtained between parametriza- 
tions of a same framework (although atypic |C<| values 
can distort the correspondence between CY\ and the fa- 
vored q values). However, it is not sufficient to explain 
the difference between relativistic and Skyrme models: 
indeed, in the relativistic case, the larger CYi are compen- 
sated by larger |C< |, leading to values of q m ax,quad similar 
to the Skyrme ones. The larger cluster sizes predicted by 
the relativistic models are due to the non-quadratic part 
of their g-dependence. 



V. STELLAR MATTER 

In the last section, we have discussed nuclear-matter 
properties in the framework of Skyrme and relativistic 
models. We now want to investigate the consequences 
of the different features we have obtained in the context 
of compact-star physics. In the first part, we discuss the 
EOS of homogeneous matter at /3-equilibrium (disregard- 
ing the liquid-gas instabilities), considering the possibil- 
ity of neutrino trapping. In the second part, we address 
the implications of the dynamic instabilities for compact- 
star properties: width of neutron-star crusts and non- 
homogeneities in the cores of type-II supernovae. 

A. Homogeneous /3-equilibrium matter 

The /3-equilibrium conditions impose the following re- 
lations between the chemical potentials of the particles: 



where fi Ue = for neutrino-free matter. Muons arc 
present if they can be in chemical equilibrium with the 
electrons, satisfying 

AV — A 4 "/* = Me — A 4 ^ ! 
the muon onset thus occurs when 

For neutrino-free matter, this condition reduces to \i e — 
/13 = m^, from which we determine the muon onset den- 
sity p^-onset for the different nuclear models. The cor- 
responding values are given in Table IV. In most cases, 
P/^-onset belongs to the interval [0.112; 0.124] fm~ 3 ; we 
can notice the higher values obtained from SGII and the 
relativistic models with 5 meson. This is a consequence 
of the symmetry-energy behavior in this range of den- 
sities. A higher symmetry energy increases the proton 
fraction at /3-equilibrium: the electron chemical poten- 
tial is then higher, and the muon onset is reached more 
easily. We can verify on Fig. 2 that, for p <~ 0.11 fm~ 3 , 
the symmetry energy curves have very similar values, ex- 
cept precisely for SGII and the models with 5 meson for 
which a s is smaller. 

In Fig. 10 top, we plot the proton fractions at (3- 
cquilibrium for neutrino-free matter. Results are shown 
taking muons into account or considering only electrons. 
As noticed earlier, the proton fraction at /3 equilibrium 



TABLE IV: Baryonic density at muon onset. For SIII the 
muons disappear at densities larger than 0.23 fm~ 3 . 



model p M -on S ct (fm 3 ) 



SIII 


0.119 


SGII 


0.146 


SLy230a 


0.121 


NRAPR 


0.117 


LNS 


0.124 


NL3 


0.112 


NL<5 


0.142 


TW 


0.115 


DDME2 


0.114 


DDH<5 


0.166 
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FIG. 10: (Color online) /3-equilibrium for homogeneous mat- 
ter at T — 0, for non-relativistic (left) and relativistic (right) 
models. Top: neutrino-free matter; the proton fraction is 
given considering the muon onset (thick line) or considering 
only electrons (thin line). Bottom: proton and neutrino frac- 
tions for a constant lepton fraction Yj = 0.4. The direct 
URCA threshold for matter without muons is indicated in 
the top figures by a thin line (DU). 



essentially reflects the symmetry-energy a s (p): see Fig. 2. 
For most models, the proton fraction increases quite 
softly with density, reaching a 1^-range of ~ [0.08; 0.11] 
at p = 0.45 fm -3 . Two kinds of atypic behaviors are ob- 
served, following the symmetry-energy features. (1) The 
two older Skyrme parametrizations show a rise and fall 
of the proton fraction with density, eventually leading to 
pure neutron matter. (2) On the opposite, the relativis- 
tic models with constant couplings show a very sharp 
increase of the proton fraction with density. DDH-<5 also 
predicts a quite sharp increase of Y p at high densities, 
due to the effect of the <5-meson on a s (p). 

However, let us consider the three modern Skyrme 
forces: SLy230a, NRAPR and LNS. Despite very close 
values of the symmetry energy, the proton fractions ob- 
tained with SLy230a are much lower. Looking at Fig. 3, 
we see this is a consequence of the different behaviors 
of the ratio p,^/ ^ axa '. For similar values of the symme- 
try energy, at high asymmetry the p% value is lower for 
SLy230a than for NRAPR and LNS: the /3-equilibrium 
is thus realised for lower Y p . The same effect can be ob- 
served among the relativistic models, comparing DDH-<5 
to the other DDH models. 

The cooling of neutron stars may occur through a di- 
rect/indirect URCA process [48]. Since the first predicts 
a too fast cooling, models which allow it are not ade- 
quate for the description of asymmetric matter. If the 



0.12 
0.1 

CO 

l c 

C 0.08 
>+— 

V) 

g 0.06 
u 

jQ 0.04 
Q_ 

0.02 






Pv-u 




SIN 




SGII 




SLy230a - 




NRAPR 


w 


LNS 




NL3 — - 



2 4 6 8 10 12 14 16 18 

T (MeV) 




8 10 12 14 16 18 

T (MeV) 



FIG. 11: (Color online) Total density at the crossing between 
the instability region and the /3-equilibrium condition, as a 
function of the temperature. Left: /3-equilibrium in neutrino- 
free matter. Right: /3-equilibrium at constant lepton fraction 
Yi = 0.4. 



muon onset is not taken into account, the critical proton 
faction that allows direct URCA is you — 1/9- In the 
presence of muons, this fraction is increased to [49] 
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Only the relativistic models with constant coupling 
present proton-fractions large enough for a direct URCA 
process in the range of density we show: the sharp a s 
evolution of these models is thus in contradiction with 
the neutron-star cooling observations. 

On the bottom of Fig 10, we consider the /3-equilibrium 
in matter with trapped neutrinos, taking a constant lep- 
ton fraction Yi = 0.4. Due to the presence of electronic 
neutrinos, for the range of density that we consider the 
muon fraction at equilibrium is vanishingly small: in this 
situation we will include only the constituents n, p, e, v e . 
On the figure, we show the proton and neutrino fractions. 
The most striking feature is that all models (except SIII) 
give similar predictions, with a nearly constant proton 
fraction in the range ~ [0.3; 0.35]: when the lepton frac- 
tion is fixed, the dependence of matter compositions on 
the symmetry energy becomes very weak. 



B. Clusterization of stellar matter 

In the present section, we discuss the formation of non- 
homogeneities in compact-star matter, considering con- 
sequences for neutron-star crust and supernova core. For 
neutron-star crust, we consider neutrino- free matter at 
T = 0; for the supernova context, we have to include 
finite temperature, and it is relevant to consider the ef- 
fect of neutrino trapping. Finite temperature results are 
given only for Skyrme models and NL3. 

Although the spinodal region almost reaches pure neu- 
tron matter at T = 0, it is limited to more symmet- 
ric matter as the temperature increases, until it disap- 
pears for a limiting value of T [40]. We wish to deter- 
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Y v = Y L = 0.4 



pcross,out pcross.in pcross,out pcross,in 

xlQ- 2 (fm- 3 ) (fm" 3 ) x!Q- 2 (fm- 3 ) (fm~ 3 ) 



SIII 


0.877 


0.110 


0.130 


0.093 


SGII 


0.418 


0.076 


0.029 


0.088 


SLv230a 


0.459 


0.079 


0.035 


0.088 


NRAPR 


0.475 


0.072 


0.031 


0.088 


LNS 


0.543 


0.077 


0.037 


0.096 


NL3 


0.553 


0.053 


0.083 


0.081 


NL5 


0.442 


0.057 


0.086 


0.090 


TW 


0.915 


0.075 


0.108 


0.084 


DD-ME2 


0.610 


0.072 


0.060 


0.083 


DDH5 


0.776 


0.079 


0.098 


0.084 



TABLE V: Predicted density at the outer (p C ross,out) and in- 
ner edge (p C ross,in) of the crust of a compact star at zero tem- 
perature, as defined by the crossing between the dynamical 
instability region and the /3-equilibrium condition for homo- 
geneous, neutrino-free stellar matter. 



mine under which conditions compact-star matter at 0- 
equilibrium reaches the region of instability against clus- 
ter formation. On Fig. 11, we plot the total density at 
the crossing between the instability region and the 0- 
equilibrium condition, as a function of the temperature. 
We call Across the maximal temperature for which this 
crossing occurs. Two cases are considered: neutrino-frcc 
matter and matter with trapped neutrinos (Yi = 0.4). 
A strong model-dependence is observed in the case of 
neutrino-free matter, where the equilibrium involves 
very neutron-rich matter. However, two common fea- 
tures can be drawn: (i) for all models, the instability 
region is crossed at T = 0; (ii) for most models, we ob- 
tain T cross < 3 MeV. With SIII, higher temperatures are 
reached because of the larger (irrealistic) instability. For 
the other parametrizations, T cross is correlated with the 
symmetry energy at low density: the highest value is 
thus obtained with SLy230a, and the lowest one with 
NL3. Indeed, higher a s values lead to a more symmetric 
composition, deeper inside the instability region. The re- 
ciprocal analysis can be made in terms of 7x3: a spinodal 
region with a large /Z3 extension is more easily reached 
at equilibrium [56] . In the case of matter with trapped 
neutrinos, a crossing is obtained until T > 10 MeV for 
all models: indeed, the proton fractions are now larger 
than 0.3. The differences we observe in the upper-density 
crossing reflect the various p s values (upper border of the 
thermodynamic spinodal for symmetric matter). 

In Table V we give the crossing densities at T = 0, 
for the neutrino-free and neutrino-trapping cases. We 
denote p C ro SS ,in (/Ocr OSS ,out) the density at the high (low) 
density crossing point. 

The value of /5 croS s.in obtained for neutrino-free matter 
at T = provides the lowest estimation for the density 
at the inner border of the crust, p C rust,in- Indeed, the 
finite-size instability region is the minimal region where 



the matter at thermodynamic equilibrium is in a cluster- 
ized shape, so it is contained by the crust. We expect 
Pcross.in to be a good approximation to p C rust,in- Most 
of the models we present give p C ross,in ~ 0.075 fm -3 . 
The lower values obtained with NL3 and NL<5 can be 
related to their larger spinodal-contour concavity C s , dis- 
cussed in Sec. IV. Nuclear matter in the crust of neutron 
stars has been studied recently [50, 51] both within RMF, 
DDH and self-consistent Skyrme Hartrce-Fock. For mat- 
ter at /3-equilibrium, the transition densities to the homo- 
geneous phase predicted in these works are: 0.085 fm -3 
(Skyrme), 0.061 fm" 3 (DDH) and 0.072 far 3 (RMF). 
These values are in reasonable agreement with the num- 
bers given in Table V, although according to our results 
with DDH models the transition density should be higher 
than ~ 0.07 fm~ 3 . We can also point out that within 
each framework the values depend on the properties of 
the chosen paramctrization. 

On the other hand, the lower-density crossing point 
Pcross.out has little significance for the crust, since the 
outer crust border does not correspond to a transition to 
homogeneous matter. At zero temperature, the very low 
density matter is always made of clusters, and the con- 
cept of homogeneous nuclear matter breaks down. Note 
however that both /9 C ross,in and /O croS s,out may be of phys- 
ical interest at finite temperature, for the formation of 
non-homogeneities in supernova cores. The /9 C ross,out val- 
ues tend to be larger with the relativistic models, which is 
another manifestation of their reduced instability region. 

Let us now consider the crossing densities for matter 
with trapped neutrinos, at Yi = 0.4. They are given at 
T = for a direct comparison with the case of neutrino- 
free matter; however, we should remark that neutrino- 
trapping occurs in the early stage of neutron-star evolu- 
tion, involving finite temperatures. The value of Across, in 
is higher than for neutrino-free matter, except for SIII: 
this also reflects spinodal-contour concavity, since with 
trapped neutrinos the equilibrium is shifted towards 
symmetric matter. The value of /9 C ross,out is typically one 
order of magnitude lower than for neutrino-free matter. 

Let us now comment on the supernova context. A 
well-known issue in type II supernova simulation is the 
difficulty to obtain the ejection of the outer layers of the 
collapsing star [52]. An additional mechanism is needed 
to produce a shock revival after the first bounce, and 
it is generally assumed that neutrino transport is a cru- 
cial factor in the explosion dynamics. It has been pro- 
posed that the liquid-gas instabilities in supernova mat- 
ter could play an important role, since it affects the trans- 
port properties of the neutrinos [53]. We are then inter- 
ested in comparing the dynamic spinodal region to the 0- 
cquilibrium condition, in the situation that can be found 
in a type-II supernova core: high temperature and pos- 
sibly neutrino trapping. In Rcf. [25], it was shown that 
an interplay occurs between neutrino trapping and clus- 
ter formation at temperatures of several MeV. Indeed, 
neutrino presence leads to more symmetric matter at 0- 
equilibrium, reaching the instability region even though 
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FIG. 12: (Color online) Comparison between the /3- 
equilibrium condition in homogeneous compact-star matter 
and the instability region in the chemical potential represen- 
tation, at finite temperature T — 10 MeV. 



it is reduced by temperature; reciprocally, cluster forma- 
tion is known to favor neutrino trapping. This feature is 
illustrated here on Fig. 12. We show the dynamic spin- 
odals at T — 10 MeV in the chemical-potential represen- 
tation (p n , p p + n e ), where the /3-equilibrium is a model- 
independent straight line, the diagonal p n = p p + p e . At 
this temperature, none of the instability regions reaches 
this line; however, all are crossed by the /3-equilibrium 
condition for a constant lepton fraction Yi = 0.4. The 
distance between each spinodal and the diagonal mea- 
sures the neutrino chemical potential which is needed to 
reach the instability region. This is related to the ps ex- 
tension of the spinodals, which gives an example of dis- 
tinction between models having high or low values of a s 
at subsaturation density: disregarding SHI, the smallest 
neutrino trapping is needed for SLy230a, and the highest 
is for NL3. 

We finally discuss the typical cluster size obtained ac- 
cording to the spinodal-instability properties, in stellar- 
matter condition. As in section IV D, we define it as the 
half-wavelenght of the most unstable mode, Ao/2. This 
quantity is given on Fig. 13 as a function of the density, 
in two different conditions: cold, neutrino-free matter 
(for all models) and hot matter with trapped neutrinos 
(for Skyrme models and NL3). For each density, we esti- 
mate the proton fraction according to the /3-equilibrium 
condition in homogeneous matter [57]. Due to the weak 
dependence of Ao on the system asymmetry, the curves 
we obtain essentially reflect the features we had obtained 
earlier for A at Y p = 0.3 (Fig. 8). At T = 0, we no- 
tice that the high asymmetry leads to a global increase 
of the cluster sizes by ~ 1 fm; a similar behavior was 
obtained in Refs. [20, 25]. The effect of temperature is 
to reduce the instability region and increase the cluster 
size (typically by ~ (2 - 3) fm for T = 10 MeV). We 
finally remark that the range of cluster sizes we obtain 



in this spinodal-scenario approach is in reasonable agree- 
ment with pasta-phases calculations performed in RMF 
framework [5, 54]. 



VI. CONCLUSIONS 

The present work provides a direct comparison of 
Skyrme and relativistic models predictions for nuclear 
and compact-star matter properties, involving the bulk 
equation of state and the finite-size liquid-gas instabili- 
ties. In our comparison, many similarities are found, and 
some differences are pointed out. 

As expected, the largest differences are obtained for the 
isovector properties. The relativistic models with con- 
stant couplings have a very hard symmetry energy. On 
the opposite, the older Skyrme forces predict a decrease 
of the symmetry energy at high density. Between these 
two extremes, modern Skyrme forces and DDH models 
present similar behaviors. However, we can notice the 
specificities of DDH<5 and SLy230a in the high density 
region: both present a stiffer a s evolution, and atypic 
rates ps/p^^ (for DDtW, this is linked to the inclusion 
of the S meson). 

Concerning the thermodynamic spinodal region, we 
have verified that its isoscalar-density extension p s is cor- 
related to the saturation density po- At zero temperature, 
the spinodal contour reaches very high asymmetry for 
all models; using the {p,ps) representation, we have ob- 
tained different p% extensions, reflecting the low density 
values of the symmetry energy. Relativistic models tend 
to yield smaller spinodal regions, both in p and p^ di- 
rections. The thermodynamic instability direction leads 
to the usual isospin distillation for all models; however 
for larger densities this effect becomes stronger for RMF 
models with constant couplings and presents a reduction 
for Skyrme and DDH models. Some very recent DBH re- 
sults seem to confirm this trend, although with a smaller 
reduction [55]. 

The dynamic finite-size instabilities leading to matter 
clusterization have been addressed in the Vlasov formal- 
ism. The wavelength associated with the most unsta- 
ble mode gives an estimation of the cluster size issu- 
ing from a spinodal decomposition: it was shown that 
Skyrme parametrizations show larger growth rates and 
favor fluctuations of shorter wavelengths. The favored 
wavelengths reflect the finite range behaviour of the force. 
In the Skyrme parametrizations, only q 2 terms are in- 
cluded. With relativistic models, the finite range is de- 
scribed by the exchange of mesons. For all models, the 
dynamic distillation effect is reduced with respect to the 
bulk one; this reduction appears stronger for relativistic 
models including the S meson. 

For the study of compact-star matter properties, we 
have considered /3-equilibrium under two different con- 
ditions: for neutrino-free matter, and for matter with 
trapped neutrinos according to a fixed lepton frac- 
tion. For neutrino-free matter, the proton fraction at 
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FIG. 13: (Color online) Half-wave-length of the most unstable mode inside the instability region. At each density, the proton 
fraction is given by the equilibrium condition for homogeneous matter. Graphs a) and b) : neutrino-free matter at T — 0, 
respectively for Skyrme and relativistic models. Graph c) : matter with trapped neutrinos (fixed lepton fraction Yi = 0.4) at 
T = 10 MeV, for Skyrme models and NL3. 



/3-equilibrium is very sensitive to the symmetry energy 
a a (p). The RMF models with constant couplings thus 
predict very high proton fractions, allowing the direct 
URCA process already at quite low densities. However, 
it was observed that the proton-fraction is not uniquely 
determined by a s (p), but also depends on the parabolic 
behavior of the isovector EOS: see SLy230a and DDH<5, 
leading to proton fractions lower than other models with 
similar symmetry energy. The situation is different for 
matter with trapped neutrinos: fixing a lepton fraction 
Yi = 0.4, all models (except SIII) predict a nearly con- 
stant proton fraction ~ 0.3 — 0.35. 

We have finally discussed clusterization of stellar mat- 
ter in two different contexts: at zero temperature, where 
it is related to the extension of the neutron-star crust, 
and at finite temperature, where it should influence su- 
pernova dynamics. In the first case, we give a lower es- 
timation (/Ocross.in) of the transition density at the inner 
edge of the crust. Our results are in reasonable agree- 
ment with values obtained within pasta-phase calcula- 
tions. Modern Skyrme parametrizations and DDH mod- 
els give similar results, p C ross,m ~ 0.75 fm . 

Stellar matter at finite temperature is addressed only 
with Skyrme models and NL3. Finite-T calculations still 
have to be performed for the relativistic models with 
density-dependent couplings, but we do not expect those 
results will affect our present conclusions. For all the 



models we show, neutrino trapping is needed to reach 
the instability region at T = 10 MeV. The required trap- 
ping rate is higher for models with low symmetry energy 
at subsaturation density, such as NL3. 

Globally, we can say that the isovector EOS shows 
quite large quantitative differences even between mod- 
ern forces: new data are still needed to better constrain 
the neutron-rich matter properties. We have shown the 
consequences of the stiff symmetry energy of the RMF 
models:stronger distillation, smaller /^-extension of the 
spinodal, and larger proton-fractions at high density, al- 
lowing the URCA process. The finite-range of the force 
is also important: its different behavior between Skyrme 
and relativistic models causes the difference in the pre- 
dictions of typical cluster size. 

A similar work comparing phenomenological models 
with results from other approaches such as Brueckner- 
Hartree-Fock should be performed. 
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